Filamin cross-linked semifiexible networks: Fragility under strain 
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The semifiexible F-actin network of the cytoskeleton is cross-hnked by a variety of proteins in- 
cluding filamin, which contain Ig-domains that unfold under applied tension. We examine a simple 
semifiexible network model cross-linked by such unfolding linkers that captures the main mechan- 
ical features of F-actin networks cross-linked by filamin proteins and show that under sufficiently 
high strain the network spontaneously self-organizes so that an appreciable fraction of the filamin 
cross-linkers are at the threshold of domain unfolding. We propose an explanation of this organiza- 
tion based on a mean-field model and suggest a qualitative experimental signature of this type of 
network reorganization under applied strain that may be observable in intracellular microrheology 
experiments of Crocker et at. 

PACS numbers: 87.16Ka, 82.35.Rs, 62.20.Dc 



The cytoskeleton of eukaryotic cells is a cross-linked 
biopolymer network^, |^ |^. Its principal constituent 
is a stiff protein aggregate (F-actin) that is cross-linked 
densely on tlie scale of its own thermal persistence length. 
Because of the combination of filament stiffness and dense 
cross-linking this semifiexible polymer gel differs funda- 
mentally from the better understood flexible polymer gels 
that are the products of modern synthetic chemistry. 

There has been considerable progress in understanding 
the complex relationship between the mechanical proper- 
ties of semifiexible networks and the mechanical proper- 
ties of their constituent filaments {4, ,5, 6, Tj & J.i IQj ■ 
Since much of this work addresses highly simplified sys- 
tems, one may ask how well does the current understand- 
ing of semifiexible networks elucidate the rheology of the 
cytoskeleton, which is a highly heterogenous chemical 
system. Cytoskeletal filaments are polydisperse in length 
and have a greater distribution of mechanical properties 
(due to e.g. filament bundling) than the model semi- 
flexible network systems studied. Furthermore these fila- 
ments are cross-linked by a plethora of highly structured 
proteins that play an active role in generating internal 
stresses and in sensing externally imposed stress. One 
class of cross-linking proteins contain numerous repeat 
domains, such as titinUl 113 and filamin[il El 

that 

unfold reversibly at a critical pulling force. 

In order to begin to address the mechanical effect of 
this chemical heterogeneity in the cytoskeleton, we inves- 
tigate networks cross-linked by filamin-like proteins con- 
taining multiple unfolding domains. We find that, above 
a certain strain threshold, the population of cross-links at 
given tension grows rapidly up to the critical unfolding 
tension of the domains. Thus, the system appears to 
self-consistently adjust its mechanical properties so as to 
reach a highly fragile state in which a large fraction of 
its cross-linkers are poised at the unbinding transition of 




FIG. 1: (color online) Model network showing the F-actin 
filaments in blue and the filamin cross-linking agents in red. 



their internal domains. The evolution of this high suscep- 
tibility state in which the system responds to small stress 
fluctuations with anomalously large strain fluctuations 
may contribute to the large nonthermal low-frequency in- 
tracellular strain fluctuations as measured by Hoffman et 
eLjl6|. The network may evolve into this high suscepti- 
bility state under the action of internal molecular motors 
(e.g. myosin - not considered in our model) so that small 
fluctuations in motor protein activity leads to the coordi- 
nated unfolding of numerous fllamin cross- linkers and the 
consequent large-scale cytoskeletal rearrangement event. 

We study via numerical simulation a random, statis- 
tically homogeneous, two-dimensional, isotropic filament 
network. These networks are formed in a manner iden- 
tical to that of Head et al.^. At filament intersections 
we add a cross-linker of zero rest length that exerts con- 
straint forces but no constraint torques. A model network 
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constructed by this procedure is shown in Fig. The fil- 
ament sections between crosslinks are modeled as linear 
springs with fixed elastic constant per unit length. The 
nonlinear behavior of semiflexible networks with freely 
rotating cross-links has been shown to be dominated by 
semiflexible filament stretching instead of bending |l7j . 
so we neglect filament bending. We anticipate that the 
results derived here are essentially independent of net- 
work dimensionality since network connectivity, not the 
dimensionality of the space in which the network is em- 
bedded, should control the collective mechanical proper- 
ties of the system. 

At forces below the unfolding force, the force extension 
relation of the filamin cross-linkers is that of a worm-like 
chain^^. When an Ig domain unfolds the contour length 
of the filamin increases, adding enough length to relax 
most of the tension at fixed extension. For simplicity, 
we model the filamin as a simple spring with spring con- 
stant kf and we take the additional contour length gen- 
erated during any unfolding event £^ to be a constant. 
The critical unfolding force of a domain is kf£f. We ne- 
glect the the rate-dependence of this unfolding force [l9|. 
Though the physiological filament cross-linkers have a fi- 
nite number of unfolding domains (24), we will take our 
sawtooth force extension curve to have an infinite number 
of branches. 

The network is sheared using Lees-Edwards boundary 
conditions. At the beginning of each strain step all nodes 
are moved affinely, then the node positions are relaxed 
through a conjugate gradient routine to a point of local 
force equilibrium. Since the cross-linker force extension 
curve is a sawtooth, there are many possible equilibrium 
states of the network. In principal, the multiplicity of 
equilibrium states requires us to use strain steps result- 
ing in displacements smaller than the sawtooth length 
^/ so that equilibrium is achieved at the smallest filamin 
extension. To reduce computational overhead we use a 
two step equilibration procedure that finds a state close 
to desired one, but allows for large strain steps. In the 
first equilibration step, we replace the sawtooth force law 
for all cross-linkers by the following force law: 

f^ffc/x |x|<^/, 
Wf |x|>^/. 

The combined network of linear elastic filaments and con- 
stant force cross-links is equilibrated. We then reimpose 
a sawtooth force law for the cross-linkers and equilibrate 
the network a second time. As the network relaxes dur- 
ing this final equilibration step, the force on each filamin 
must be less than kf£f, so the cross-links will stay on the 
same sawtooth branch. Since the rest of the network was 
originally equilibrated at the critical pulling force, the 
sawtooth force law could not have reached force equilib- 
rium on any earlier sawtooth branch assuming all filamin 
linkers act independently. In practice, collective relax- 
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FIG. 2: (color online) Distribution of normalized cross-linker 
lengths £/£f modulo 1 in equilibrated networks with, from 
shallowest to steepest slopes respectively, kf — 100, 200, 600, 
2000, 4000, and 10000. 

ations of the network push individual cross-links onto 
different sawtooth branches in this final step. We found, 
however, that such coordinated relaxation events had a 
negligible quantitative effect on the data. 

We present data for networks composed of monodis- 
perse filaments of length in = 0.2 (all lengths are mea- 
sured in units of the length of the square simulation box) 
at a filament density such that there are on average 30 
cross-links per filament. For these values we find negli- 
gible system-size effects. The length of the filamin do- 
mains is given by £f = 1.3 x 10^^, or £//4 = 0.02. This 
ratio is about ten times smaller than the expected phys- 
iological value 1^ . We choose the smaller value of 
£f/£c because it enhances the effects we were measur- 
ing; our initial studies at physiological values of £f/£c 
find qualitatively similar results, but the onset of non- 
linear effects occurred at higher strain values To 
fix an energy scale we set the extensional modulus fi of 
the filaments to unity. The average spring constant for a 
filament segment can then be determined from the mean 
distance between cross-links: ka = l/(^c) = 150. The 
range of cross-linker spring constant values studied here 
is 10^ <kf < 10". 

Figure |21 shows the measured equilibrium distribu- 
tions of cross-link lengths, modulo the sawtooth length 
£f, for a representative set of strained networks with 
£f = 1.3 X 10"'' and several values of spring constant 
kf. For values oi kf < 10 x (ku) the statistical weight 
for finding a cross- link extension (modulo £f) is expo- 
nentially enhanced towards length £ f where the domains 
unbind. For values kf > 10 x (kn) the statistical weight 
for finding a cross- link extension (modulo £f) grows faster 
than exponentially near length if. The enhancement of 
the number of such filamins poised at the unbinding tran- 
sition (critical cross-linkers) is one of principal results of 
our work. Apart from the noise floor, the shape of the 
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distribution appears to be strain independent. 

To understand the growth of the critical cross-Hnker 
population, we consider a mean-field model for the be- 
havior of a single cross-link in an effective elastic medium 
representing the rest of the network. The surrounding 
effective medium is modeled as a single harmonic spring 
with spring constant k. Reflecting the network structure, 
the cross-linker is connected in series with the effective 
network spring. We set the total strain on the two springs 
in series (by fixing their total length) so that the cross- 
link has crossed at least one branch of the sawtooth func- 
tion. Upon the further application of extensional strain, 
the two springs in series will both linearly increase their 
extension until the filamin spring with spring constant 
kf reaches is maximum force kf£f where it is poised at 
the top of its saw-tooth force-extension curve. 

Now consider an infinitesimal increase in the total ex- 
tension that drives the unfolding of one more filamin do- 
main. Before the extension the two springs were in force 
balance so that kfif = kx where x represents the ex- 
tension of the medium spring. After the infinitesimal 
extension, the system achieves force balance on the next 
branch of the saw-tooth filamin force-extension curve so 
that extension of the filamin spring is now increased by 
if — d while the extension of the medium spring is de- 
creased to X — {if — d). Force balance requires that d, 
the distance between the current extension of the filamin 
spring and the edge of the next saw-tooth, is given by 
d{k) = kif / {k + kf). In other words, the combined sys- 
tem once equilibrated with the filamin spring at its max- 
imal force is now equilibrated with that filamin spring on 
its next saw-tooth branch at a smaller force. The strain 
in the surrounding medium has also decreased due to the 
extension of one more filamin domain. 

To maintain force balance, the filamin spring cannot 
relax its length more than if — d. Upon further exten- 
sion the filamin spring will only extend until another do- 
main unbinds. Thus in steady-state the filamin spring 
will evenly sample all extensions (modulo if) between 
if — d{k) and if. For a given value of the spring constant 
of the medium we expect that the extensions (modulo if) 
of the filamin cross-linkers Xf to be uniformly distributed 
between the bounds given above so that this distribution 
can be written as 



P{Xf,k) 



1 



d{k) 



e {xf ~ [if - d{k)]) , 



(2) 



where is a step-function. Different cross-links in the 
network, however, will not have the same local environ- 
ments; the values of k will be sampled from some statisti- 
cal distribution K (k). Integrating over that distribution 
we write the probability of finding a given filamin length 
(modulo if) Xf. 



P{Xf) 



k + k 



Ifk 



^K{k)dk. 




FIG. 3: Distribution of local effective spring constants, sam- 
pled on small sets of highly stretched filamin cross- linkers with 
kf = 600. The solid line is a fit to Eq.|I|with k/k set to 7.3. 



The step function fixes the lower limit on the /c-integral. 

We examine the distribution of the local spring con- 
stants in the random network and concentrate on the 
high-fc tail of that distribution. One may imagine that 
the effective spring constant representing the medium can 
be represented as a small number of chains of springs. 
Each chain of springs represents one path for force prop- 
agation through the random network; it is made up of 
a large number of statistically independent springs con- 
nected in series. In order to find an extremely large value 
of the effective spring constant k it must be that for one 
of the force paths all of the constituent spring constants 
are large, since the compliance of the springs in series will 
be dominated by any single soft spring. We expect the 
probability of such a rare event to be Poisson distributed 
so that, in the high-fc tail, the distribution K{k) takes 
the form 



K{k) - (fc)e-'^'/'= 



(4) 



where H is some regular function characterizing the 
small-fc behavior of the distribution {H{x) — *■ const as 
X — > oo) and the constant fc is undetermined. Combining 
Eas. I3l4l we find that P{xf) takes the form 



k / kf{xf — 

if V ^Xf 



if V ' kxf 



(5) 



as long as fc 



f^-^-^ is large enough that K{k) within the 
integral in Eq.|3can be replaced by its high-fc asymptotic 
form. Eq|5l shows the sought after exponential peak as 

Using Eq. El we may determine the ratio fc//fc using 
(3) the slope of the numerically measured distribution oi Xf 
shown in Figure |21 From fitting the data for fc/ = 600 we 
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find that kf /k = 7.3. By numerically sampling the local 
mechanical response in many realizations of the network, 
we independently verify the principal physical insight in 
the above discussion: for small values of fcj , K{k) appears 
to have an exponential tail in the stiff spring limit. This 
data is presented in Fig. El for kj — 600. The fitted hne 
demonstrates that the high k tail of this data is consitent 
with a value of fc/fc = 7.3. The agreement between the 
two independently determined values of k supports our 
analysis. In upcoming work, we demonstrate that for 
higher kf the distribution function K{k) takes a different 
form, but it is still related to the distribution function 
P{xf) via Eq. 131 "251. This latter equation should apply 
for general distributions K{k). 

We have found that the strained filamin cross-linked 
network develops into a highly fragile mechanical state in 
which a large fraction of the cross-linking filamins reach a 
critical state where they are poised at the brink of domain 
unfolding. The presence of fluctuating internal stresses 
in the cytoskeleton produced by molecular motor activ- 
ity and/or equilibrium fluctuations can act on this highly 
fragile state to produce large strain fluctuations due to 
the correlated failure of numerous critical cross- linkers. 
Thus, the observation of the formation of this critical 
state under applied stress may explain a particular fea- 
ture of the observed low-frequency strain fluctuations as 
observed by intracellular microrheology. 

We have presented a simple, mean-field theory to ex- 
plain the evolution of this fragile state under applied 
strain. There are a number of extensions of this work 
that remain to be considered. Foremost among these is 
the exploration of the effect of filamin-type cross-linkers 
in semiflexible gels where the filaments each have a fi- 
nite bending modulus. In addition, the development of 
a complete model that includes the effect of internally 
generated random stresses due to the action of molecu- 
lar motors will be an important step towards the direct 
calculation of the low- frequency dynamics of this biopoly- 
mer gel of fundamental biological importance. 
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